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ON  SLOW  VISCOUS  FLOW 
by 

George  F.  Carrier,  Brown  University 


ci 

1,  Introduction.  The  flow  of  a viscous  fluid  past  an 

obstacle  at  low  Reynolds  Number  has  been  the  subject  of  numerous 

investigations.  A rigorous  analysis  of  the  motion  of  the  fluid 

requires  the  solution  of  a non-linear  problem  but  various  attempts 

have  been  made  to  formulate  an  equivalent^"  linear  problem.  The 

2 

work  presented  here  is  a discussion  of  a new  linearization  of 
this  problem  which  is  based  on  a conjecture  rather  than  on  any 
formal  procedure  such  as  a perturbation  process.  The  analyses 
of  several  specific  boundary  value  problems  using  this  lineariza- 
tion are  presented  and  the  results  are  compared  to  those  of  the 
clasr ical  theories  of  Stokes  [l]  and  Oseen  [2].  The  results  re- 
ported here  are  in  better  agreement  with  the  physical  facts  than 
those  given  by  the  Stokes  or  Oseen  theories.  We  shall  compare 
the  point  of  view  adopted  in  these  classical  methods  with  that 
adopted  in  ours  as  the  analysis  proceeds. 

The  actual  details  of  the  various  problems  treated  here 
were  carried  out  by  several  people  at  Brown  University;  the  ex- 
plicit references  are  included  in  the  appropriate  section.  How- 
ever the  author  is  especially  indebted  to  Professor  G.  W,  Morgan 


1.  "Equivalent"  in  the  sense  of  yielding  a good  approximation 
to  the  rigorous  solution. 

2.  This  linearization  was  actually  used  earlier  on  a specific 
problem  by  J.  A.  Lewis  and  the  author  [3]. 
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who  provided  both  the  material  and  presentation  of  section  6. 

2.  The  slow  viscous  flow  problem*  The  motion  of  a 
viscous  incompressible  fluid  is  governed  by  the  familiar  laws 
requiring  the  conservation  of  momentum  and  mass.  For  an  isotropic 
homogeneous  medium  of  constant  viscosity  ( i. e. , p.  independent  of 
thermodynamic  state)  these  laws  take  the  form 


UJUi,J  + p"  p,i 


= vu 


i,jj 


(2.1) 


UJ,J  = °‘ 


(2.2) 


Here,  u.p  p,  p,  v,  are  respectively  the  velocity,  pres- 
sure, density,  and  kinematic  viscosity.  The  differentiations 
are  performed  with  regard  to  the  physical  coordinates  x{. 

The  boundary  conditions  which  are  typical  of  some  of 
the  problems  of  physical  interest  require  that  the  velocities  bo 
specified  on  the  boundary  curves  of  the  region  under  considera- 
tion. The  type  of  problem  associated  with  such  conditions  is 
that  in  which  we  are  interested  here. 

The  differential  equations  can  be  put  in  a dimension- 
less form  by  introducing  the  following  substitutions:  Wj_  = u^/uQ, 
uQ/ v,  o = p/pu^.  Here,  uQ  is  a characteristic  velocity 
of  the  problem,  e.g, , the  free  stream  velocity.  Equations  (2.1) 
and  (2,2)  then  take  the  form 


Wi,j  + V = wi,jj 
WJ,J  = 0 


(2.3) 


(2.4) 
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where  the  differentiations  are  now  performed  with  regard  to  the 
x^.  This  is  the  form  of  the  equations  which  we  shall  find  con- 
venient but,  in  order  that  we  may  compare  the  classical  analyses 
of  this  problem,  we  note  that  another  formulation  is  readily  ob- 
tained when  one  else  introduces  a length  a which  is  characteris- 
tic of  the  geometry  of  the  problem*  We  can  then  define  e = uQa/v , 
= x.|/a,  v^  = u^uQ,  and  x = p/pUg.  Then  Eqs.  (2.1)  and  (2.2) 


take  the  form 


evj  vi,j  + V = vi,jj 


= ° 


{2.5) 

(2.6) 


where  the  differentiations  are  taken  with  respect  to  the 

Stokes,  who  was  concerned  with  the  flow  of  an  otherwise 
undisturbed  uniform  stream  past  a solid  motionless  sphere,  intro- 
duced a perturbation  procedure  using  e (the  Reynolds  number)  as 
the  expansion  parameter.  That  is,  he  used  the  representations 

= #0  w*"  * = g 


and  obtained  a sequence  of  linear  problems  for  the  determination 
of  the  Vjn  and  xn  when  these  series  were  introduced  into  Eqs. 
(2.5)  and  (2.6)  and  into  the  boundary  conditions.  When  e is 
sufficiently  small,  the  results  of  such  an  analysis  are  in  excel- 
lent agreement  with  the  experimental  observations  (Fig.  1)  for 
the  flow  past  a sphere.  However,  when  this  procedure  is  attempt- 
ed for  certain  two  dimensional  problems,  it  is  found  that  no 
solutions  of  the  linear  problems  so  formed  can  exist.  Specifi- 
cally, those  problems  in  which  we  require  that  u^  = 0 on  a closed 
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curve  and  u^— > uQ,  U£ — >0  as  the  distance  from  this  curve  in- 
creases do  net  admit  a perturbation  series  of  this  type.  This 
fact  led  Oseen  to  modify  the  Stokes'  procedure  somewhat.  Essen- 
tially, using  the  notation  of  Eqs.  (2.3)  and  (2,4),  he  developed 
the  velocity  components  as 

= 1 + w^_  + 

I 

W2  - W2  "f  * , , , 

-J 

where  the  primed  quantities  are  of  "higher  order"  . That  is,  he 
wrote  the  velocity  field  as  the  undisturbed  velocity  plus  the 
perturbation  induced  by  the  obstacle.  In  each  of  these  proce- 
dures, one  computes  only  the  terms  of  lowest  order  and  these 
give  results  whose  macroscopic  features  are  reasonably  accurate 
approYimatlons  to  the  rigorous  result  for  e £ 1, 

An  alternative  way  of  viewing  these  foregoing  proce- 
dures (where  only  the  first  term  of  the  expansion  is  found)  does 
not  require  the  formal  introduction  of  a perturbation  process  or 
successive  approximation  scheme.  One  may  consider  an  approxima- 
tion technique  where,  in  Eq.  (2.1),  we  replace  Uj  by  sero  (for 
the  Stokes  model)  and  u0  6^  (for  the  Oseen  model).  In  the  nota- 
tion of  Eq.  (2.3)  then,  we  have 

°3  wi,3  + °,i  = wi,jJ  (2*7) 

3.  The  detailed (argument  can  be  found  in  [2],  It  is  clear,  how- 
ever, that  wi  must  take  on  the  value  -1  on  the  obstacle  and 
hence  that  w7  is  not  strictly  of  higher  order.  Thus,  Oseen* s 
formulation  can  be  taken  as  the  basis  of  a successive  approx- 
imation scheme  but  not  as  a formal  perturbation  procedure 
in,  say,  e . 
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where  Cj  = (0,0)  or  (1,0)  for  the  two  models  in  question. 

In  order  to  motivate  the  forthcoming  suggestion  for  a 
modification  of  these  linearizations  it  is  convenient  (out  not 
necessary)  to  consider  the  two  dimensional  problem.  When  we  re- 
strict ourselves  in  this  way,  we  can  write  w-^  = i|jy,  w~  = “ 
and  Eo.  (2.3)  implies  that 

= i|>y  - 4>x  A4y  (2.8) 

When  we  use  Ea.  (2.7)  instead  of  (2.3)  we  obtain 

= cx  A\|»x.  (2.9) 

Let  us  now  consider  obstacles  with  a sharp  leading  edge 
which  is  positioned  at  the  origin.  The  stream  function  i|>  which 
is  to  satisfy  the  boundary  conditions  and  Eqs.  (2.8)  or  (2.9) 
(depending  on  whether  we  choose  the  exact  problem  or  the  approxi- 
mate one)  must  have  a branch  point  at  the  origin.  One  sees  this 
in  the  following  way.  Let  the  boundary  curve  y(x)  be  an  analytic 
function  near  and  at  the  sharp  leading  edge  and  consider  the 
function  [x,y(x)3  on  this  curve  and  on  its  analytic  continua- 
tion into  the  flow  field.  Then  on  the  boundary,  u^  [x,y(x) ] is 
identically  zero  but  on  the  continuation  it  is  not.  It  follows 
that  is  not  analytic  at  the  origin.  That  is  must  behave 
like  rn  near14-  r = 0 and  n must  lie  between  the  values  3/2  and 
2,  [r  = (x^  + x|)*].  For  a cuspital  edge  (see  [3])  n is  3/2. 

This  implies  that  the  contributions  of  the  fourth  derivatives  on 
the  left  of  2q.  (2.8)  [or  (2.9)]  are  more  strongly  singular  than 
those  terms  on  the  ric.ht*  This  further  implies  that  is  a func- 

tion  whose  behavior  near  r = 0 is  predominantly  that  of  a 

The  singularity  could  be  worse  than  this,  depending  on  the 
leading  edge  shape. 
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biharmonic  function.  That  is  to  say,  if  we  subtract  from  this 
singular  biharmonic  function  the  biharmonic  term  will  be  more 
highly  singular  than  ^ - D and,  in  fact,  -4»/D  — > 1 as  r — >0.  Here 
D is  the  biharmonic  function  of  order  r11  Just  mentioned. 

> Similarly,  we  can  argue  that  as  r— > oo  , the  viscous 

terms  are  not  of  importance  and  should  be  like  a non-vi scous 
flow.  That  is,  \|>  should  become  harmonic  as  r — Note  that 
either  Eq.  (2.8)  or  (2.9)  admits  such  harmonic  solutions  (any- 
where). The  important  fact  to  note  now  is  that  while  each  equa- 
tion together  with  the  boundary  conditions  is  capable  of 
"generating11  a function  ^ which  is  predominantly  blharmonic  near 
. r = 0 and  harmcnlc  as  r— »co  . their  behaviors  in  the  intermediate 

region  are  not  identical  and  there  is  little  reason  to  hope  that 

% 

each  equation  will  imply  a ^ which  continues  from  ( sav)  a given 
blharmonic  behavior  near  r = 0 to  the  same  harmonic  function  at 
large  r.  However,  it  is  to  be  hoped  that  if  were  neither  1 
nor  0T  but  had  a value  corresponding  to  some  appropriately  weight- 
ed average  of  the  velocity,  the  differences  in  the  values  of  the 
ao-lat-ions  of  the  two  equations  (2.8)  and  (2.9)  would  vanish  on 
the  average.  If  this  should  occur,  then  Eq.  (2.9)  and  its  bound- 
ary conditions  could  generate  a solution  wnich  was  essentially 
correct  at  large  r,and  near  the  front  of  the  obstacle,  but  which 
was  not  necessarily  a very  accurate  approximation  in  between. 

Since  the  boundary  conditions  are  to  be  applied  where 
the  solutions  are  hypothetically  correct,  however,  the  macro- 
scopic features  of  the  prediction  would  be  acceptable. 

We  should  note  that  if  this  sort  of  an  "average 

i s 
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equation",  (2,9),  should  be  acceptable  as  a substitute  for  (2,8), 
it  is  probable  that  c^  would  be  a function  of  the  important  dimen- 
sionless parameter  appearing  in  such  problems,  i,e. , the  Reynolds 
number  e. 

Our  conjecture  then  is  that  there  may  exist  a function 
c(e),  0 < c < 1,  such  that  the  solution  of  Eq.  (2,8)  under  suit- 
able boundary  conditions  can  be  "replaced"  by  a solution  of  Eq. 
(2.9)  with  ci  replaced  by  c(e).  If  the  conjecture  is  to  be  use- 
ful it  is  necessary,  of  course,  that  c depend  only  on  e. 

In  the  succeeding  sections  we  shall  present  the  manner 
in  which  we  chose  the  presumably  correct  value  of  c for  moderately 
small  e and  the  comparison  of  the  results  of  this  theory  with 
either  experimental  observation  or  the  solution  of  the  corres- 
ponding non-linear  problem.  It  is  to  be  noted  that  although  we 
have  based  our  conjecture  on  a motivation  associated  with  a sharp 
edged  obstacle,  we  shall  "test"  it  on  a broader  group  of  problems. 

3«  The  flow  oast  a long  .^lat  plate.  The  classical 
problem  associated  with  the  laminar  flow  of  a viscous  fluid  past 
a semi-infinite  flat  plate  is  the  problem  which  gave  rise  to  the 
investigations  presented  in  this  paper.  In  1948,  C.  C.  Lin  and 
the  author  attempted  a calculation  of  the  flow  in  the  neighborhood 
of  the  leading  edge  of  such  a plate  [4],  The  critical  argument 
of  that  investigation  required  that  the  region  in  which  the  lead- 
ing edge  solution  was  valid  overlap  the  region  in  which  the 
Blasius  solution  was  valid.  However,  there  seems  to  be  no  con- 
vincing argument  that  such  an  overlap  region  exists.  Before  the 
doubtful  validity  of  this  result  was  noticed,  however,  the  same 
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problem  was  considered  [3 3 using  the  Oseen  formulation  and  at 
that  time  the  conjecture  of  this  paper  was  introduced.  In  order 
to  force  the  result  of  this  linearized  problem  to  match  the  lead- 
ing-edge result,  the  number  c was  given  the  value  .35*  The  solu- 
tion of  the  linearized  problem  is  outlined  in  the  appendix  (the 
details  can  be  found  in  [31)  and  the  result  is  given  by 

2 

*P(x,y)  = 2£[q  erf  (qc1//2)  - (me  )“1//2(1  - e“Cy  )] 

= S g(q)  (3-D 

1/2 

where  £ + iq  = (x  + iy)  and  tp  is  the  stream  functions  defining 
the  flow. 

When  c = .35,  the  velocity  gradient  at  the  plate  is  the 
same  as  that  found  by  Blasius  but  the  behavior  for  l^rge  q (and 
any  O differs  from  that  of  the  Blasius  solution.  In  fact,  if 
we  denote  the  Blasius  result  by  tPg  and  ours  by  <P  then 

?B  = 5f(q)  (3.2) 

and  asymptotically  (for  large  q) 

9>b^  t,  (2q  - 1,72  + ...  ) 

V ^ K [2q  - 2Uc)~1/2  + ...  ] (3.3) 

These  asymptotic  developments  agree  only  when  c = .43. 

It  would  now  appear  that  one  might  choose  c = ,35  if 
the  Lin-Carrier  result  were  believod  to  be  valid,  or  c = .4-3  if 
it  seemed  appropriate  to  match  the  large  q behavior  of  the  linear- 
ized answer  and  the  Blasius  result.  The  latter  can  be  argued 
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against  by  noting  that  the  Blasius  result  is  the  first  term  of  an 
asymptotic  development  in  £ and  hence  is  valid  only  for  £ » 1 
whereas  one  can  anticipate  the  validity  of  our  linearized  result 
only  for  small  The  argument  against  the  former  has  already 
been  noted. 

Before  the  foregoing  discrepancy  was  noted,  we  proceed- 
ed to  investigate  the  problems  presented  in  sections  (4),  (5), 
and  (6).  For  moderataly  low  values  of  the  Reynolds  number,  the 
results  for  the  problems  of  these  sections  provide  excellent 
agreement  between  the  non-linear  theory  or  experimental  evidence 
and  the  c = .43  result.  This  fact,  however,  does  not  justify 
using  .43  (or  any  other  c)  at  *dgh  Reynolds  number. 

In  order  to  determine  whether  any  value  of  c might  be 
appropriate  for  the  high  Reynolds  number  problem  of  this  section, 
the  flow  field  for  the  semi-infinite  plate  problem  is  being  com- 
puted from  the  non-linear  formulation,  via  the  relaxation  method. 
The  investigation  has  not  yet  been  completed  but  it  is  clear  that 
the  local  linearised  prediction  for  c = .43  will  not  agree  with 
the  non-linear  prediction.  This  leads  to  the  amusing  situation 
where  the  choice  c = .43  which  is  very  good  at  small  and  moderate 
values  of  the  Reynolds  number  resulted  from  the  analysis  of  an 
infinite  Reynolds  number  problem  to  which  it  is  not  particularly 
appropriate. 

4.  The  .flow  past  the  cylinder  and  the  sphere.  The 
application  of  the  present  theory  to  the  determination  of  the 
flow  of  an  otherwise  undisturbed  uniform  stream  past  an  obstacle 
is  a rather  simple  matter  when  the  corresponding  Oseen  result  is 
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available.  A discussion  of  the  flow  past  a sphere  will  illus- 
trate this  fact  clearly  and  the  extension  to  other  obstacles 
(including,  of  course,  the  cylinder)  is  a straightforward  matter. 
In  this  discussion  it  is  most  convenient  to  use  the  notation  of 
Eq.  (2. 5) • In  our  linearized  theory  the  momentum  equation  takes 


the  form 


cev.  , + a . 

,1 


(4.1) 


The  Oseen  theory  differs  from  this  in  that  c is  re- 
placed by  1.  The  boundary  value  problem  requires  the  solution 

of  Eqs.  (2.6)  and  (4.1)  subject  to  the  conditions  vj  = 0 on 
2 

r = XjXj  = I and  Vj — as  r — ^°°  • 

We  now  denote  the  tangential  velocity  component  asso- 
ciated with  the  solution  of  the  Oseen  problem  (c  =1)  by  v(e,r,©) 
and  that  associated  with  the  modified  linearization  as  V(e,r,©). 
It  is  readily  seen  that  V(e,r,©)  = v(ce,r,©).  The  other  velocity 
components  and  the  pressure  may  be  related  in  the  same  way*  The 
friction  drag  on  the  sphere  (according  to  the  present  theory)  is 


given  by 


0 = 2 up.  uQa 


V (e ,1,©)  sin  © d© 


<JG 

u 

n 


= 2 xp  u-a 


v (ce,l,©)  sin2©  d© 


0 


and  the  corresponding  drag  coefficient  is  given  by 


CD  = (p/puQa)  vr(ce,l,©)  sin2©  d©. 


However,  the  friction  drag  coefficient  C^,  given  by  the  Oseen 
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theory  for  Reynolds  number  ce  is 

it 

"i 

CD  = « s1"29  d« 

D 

so  that 

CD  = c C^Cce).  (4,2) 

The  contributions  to  the  drag  of  the  pressure  variation  over  the 
sphere  are  related  in  the  same  './ay.  Thus,  Eq,  (4,2)  expresses 
the  drag  coefficient  relationship  for  the  two  theories  for  the 
gross  drag  or  for  either  contribution.  The  same  formula  is  valid 
for  the  cylinder  or  any  other  obstacle  family  characterized  by  a 
single  length  parameter. 

In  view  of  this  result  we  may  use  the  results  of  drag 
coefficient  calculations  of  previous  authors  to  find  those  asso- 
ciated with  our  theory.  Figure  (4,1)  indicates  the  results  of 
this  calculation  for  the  cylinder  using  the  figures  obtained  by 
Eairstow  [8]  and  also  indicate  some  experimental  evidence  con- 
cerning the  drag.  Figure  (4,2)  indicates  the  corresponding  evi- 
dence for  the  sphere.  It  is  clear  that  the  present  theory  with 
c = ,43  gives  a much  more  accurate  prediction  of  the  drag  than 
do  previous  accounts  of  the  matter. 

5.  The  finite  flat  Plato,  In  this  section  the  low 
Reynolds  number  flow  past  a flat  plate  will  bo  discussed.  The 
formulation  of  the  problem  differs  from  that  of  the  appendix  only 
in  the  respect  that  the  boundary  condition  = - 1,  on  y = 0, 
x > 0,  applies  only  when  x < a.  Here  a is  the  length  of  the 
plate  in  units  v/u0« 
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The  analysis  of  this  problem  follows  that  of  the  appen- 
dix up  to  Bq.  (3.4).  At  that  stage  of  the  preceedings,  it  is 
convenient  to  note  that 

7(0 


y 2(C-lk)1/2[(^+ik)1/2+({;+i)1/2] 


= f (Oh(^)  (5.1) 


and  to  apply  the  inversion  integral  to  each  side  of  this  equation 
using  the  convolution  theorem  to  evaluate  the  right  hand  side. 

We  obtain^  a 

P 

f(x)h(x  - *t)dx 


where 


'I’ytX.O)  = 


h(x)  = (£) 
XX 


cx, 


i2 


eu  kJ  uj  du. 


0 


Since  f(x)  must  have  a singularity  of  order  I at  x = 0 and  at 
x = a,  it  is  convenient  to  write 


f(x)  = g(x)/V/x(a  - x) 

* 


4»y(x,o)  = 


g(x)[h(x  - x )/  V/r  (a  - t )]dr. 


0 


2 2 

The  substitutions  t = a sin  9,  x = a sin  9,  lead  to  the  equation 


x/2 

n 


^y(x,0)  = 1 = 


2 

g(a  sin  9)  h (a  sin29  - a sin29)d9* 


This  integral  equation  for  g may  now  be  solved  either 
by  finding  the  coefficient  of  the  expansion 

g(a  sin29)  = S an  9n 


T*  A more  detailed  account  of  this  analysis  can  be  found  in  [6]. 
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or  by  a numerical  process.  It  turns  out  that  an  excellent  approx- 
imation to  g for  moderate  values  of  a is  given  by 

g(a  sin2©)  = A(a)  + B(a)(©  - u/4)« 

The  quantities  A and  B are  plotted  in  Fig.  (5»1)»  In  particular, 
the  drag  is  given  by  p,U07iA/2.  The  flow  past  a finite  flat  plate 
according  to  the  non-linear  theory,  Eqs.  (2.1)  and  (2.2),  has 
been  found  by  Munier  [7l*  Using  the  relaxation  technique  the 
flow  field  [i. e.  , Uj_(x,y),  p(x,y)]  has  been  computed  for  a = 4 
and,  in  particular,  the  ’-elocity  gradient  at  the  plate  has  been 
recorded.  This  velocity  gradient  as  deduced  from  each  of  these 
theories  is  given  J.n  Fig.  (5*2).  The  agreement  of  the  results 
of  the  linearised  theory  with  those  of  the  non-linear  theory  is 
surprisingly  good, 

6,  The  flow  in  a v/edee  shaped  region.  In  this  section 
we  shall  present  the  linearised  analysis  of  the  converging  flow 
in  a wedge  shaped  region  and  compare  the  results  to  the  exact 
solution  obtained  by  Hamel.  This  problem  is  not  in  the  category 
for  which  the  linearization  is  designed.  However,  it  is  of  inter- 
est to  see  how  well  the  flow  is  predicted  by  this  theory. 

We  start  from  Eqs.  (2.1)  and  (2.2)  and  use  the  radial 
coordinates  r,©.  The  flow  proceeds  towards  the  origin  in  the 
wedge  shaped  region  - ©0  < © < ©0»  We  look  for  solutions 
iji  = ^(©)  corresponding  to  which  the  radial  velocity  is  given  by 
u(r ,9)  = r”^\j)’(©)  and  the  circumferential  velocity  vanishes 
everywhere.  We  denote  by  e (the  Reynolds  number)  the  quantity 
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URp/p.  where  u(r,9)  = 


- U -S  and  u(r,9)  can  be  written 
r 


u(r ,9)  = - (UR/r)x‘(9), 


Equation  (2.1)  and  (2.2)  new  yield 


AAX  + 6r"1x,(  ^x)r  = 0. 


(6.1) 


Since  x is  a function  of  9 only 


X + 4x"  - 2ex,xn  = 0. 


(6.2) 


Note,  however,  that  in  Eq.  (6.1),  the  quantity  xf/r  corresponds 
to  tne  coefficient  Vj  in  Eq.  (2.5)  and  hence  corresponds  roughly 
to  the  terra  which  is  replaced  by  unity  in  the  Oseen  treatment 
and  by  c in  ours.  However,  it  is  profitable  to  retain  the  infor- 
mation that  u behaves  like  1/r  and  to  replace  by  c retaining 
the  1/r  contribution.  This,  clearly,  is  in  the  spirit  of  our 
approximation  technique  since  x'(®)  varies  from  unity  to  zero  and 
has  some  "average",  c.  Our  linearized  equation  then  takes  the 


form 


X + ^X"  ~ 2cex"  = 0. 


(6.3) 


6.1. 


■t.  Equation  (6.2)  can  be  inte- 


grated twice  to  yield 


+ 2w2  - i w^  + Cw  = D 

3 


(6.4) 


where  C and  D are  arbitrary  constants  which  must  be  determined 


and  we  have  put  w = x’»  We  have  x*(0)  = 1 and  therefore  w(0)  = 1. 
Also,  because  of  the  symmetry  about  9 = 0,  we  must  have  w'(0)  = 0. 


These  two  conditions  together  with  Eq.  (6.4)  give  the  following 


relation 
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2 - 1 C = D.  (6 .5) 

Hence  (6.4)  becomes 

(w*)2  +2w2-iw3+Cw=C--|  + 2*  (6.6) 

2 3 3 

To  integrate  Eq.  (6.6)  we  rewrite  it  as  follows: 

2 

ly.M-  = i.  (w  - l)[w2  + (1  - £)v  ->  (1  - + 6)]#  (6.7) 

2 3 e e e 

We  now  factor  the  expression  in  square  brackets  so  as  to  obtain 
(6.7)  in  the  form 

(v1)  = | e(w  - l)(w  - a1)(w  - a2)  (6.8) 

where 

ax  = \ f( | - 1)  + l/(l  - |)2  + It  (30  + 6 - e)J 

a2  = H(t  ' X)  " ^ Q - f>2  (30  + 6 - e)j 

Using  Eq.  (6.6)  and  the  boundary  condition  of  viscous 
flow  that  w(  + ©Q)  = 0,  we  see  that  C - -|  + 2 > 0 and  hence  that 

(X  - 4>2  + i Oe  + 6 - e)  > o. 

Thus  a^  and  a2  are  real. 

Equation  (6.8)  can  now  be  integrated  in  terms  of  ellip- 
tic functions.  We  have  from  (6,8) 

w*  = + yfl'  e V(w  - l)(w  - a1)(w  - a2)  . (6.9) 

Since  w(0)  = 1 and  w(  + ©0)  =0,  we  expect  that  w will  range  between 

0 and  1 and  that  w'  J>  0 between  -w0  and  0 and  w'  _<  0 between  0 

and  ©c.  Hence  the  positive  sign  in  (6,9)  holds  for  © _<  0. 
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Let  us  Integrate  from  -©0  to  some  non-positive  angle 


©. 


Va  - 1)U  - ax)U  - a2)  JQo 


dn  = o + e.  (6.io) 


© < o 


Now  a^  and  a2  still  contain  an  unknown  constant  C which  must  be 

evaluated  from  the  condition  that  w = 1 and  © = 0.  Hence 

1 

r> 


£ 


d* 


V(5  - DU  - a1)(?  - a2) 


= ®o* 


(6.11) 


Upon  integration,  Eq.  (6.11)  will  furnish,  in  terms  of  elliptic 
functions,  a formula  relating  C,  e,  and  ©Q.  Substituting  this 
into  (6.10)  and  integrating  we  obtain  the  solution  of  the  problem. 

Tables  of  elliptic  integrals  give  the  following  relation 

d£  ? -1 


R 


w 


V(o  - 5)(P  - <)«  - r)  \J  a . r 


sn 


a~Y  P-w  . /p-Y 
P-Y  a-w  V a-Y 


(6. 12) 


a > p > w > y • The  inequality  is  satisfied,  as  we  shall  soon 
see,  if  we  let  a-^  = a,  1 = P,  &2  - Y.  Now 


1 

n 


•J  -#»•- 

w 


* 


w 


by  using  Eq,  (6,11),  Hence,  using  (6.10),  we  have: 


dS 


t V^ai  " 0(1  - OU  - a2} 


V ©. 


(6.13) 
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We  now  consider  three  special  cases  corresponding  to 
different  ranges  of  the  Reynolds  number,  e • 

(a)  e < 6;  special  case  e = 1. 

For  e = 1,  the  roots  a^,  a^  are 

^ {5+^5  + 12c  } 

a2  = \ + 12o  } 


We  previously  saw  that  C > - 2.  Hence  a^  > 5*  a2  < 0,  and  the 

inequality  in  (6  12)  is  satisfied. 

(b)  e = 6. 

In  this  case  the  roots  are 


al 


To  estimate  the  values  of  a-^  and  a2  we  make  use  of  another  in- 
equality containing  C.  Since  u(r,0)  will  have  a minumum  at 
© = 0,  w(0)  will  have  a maximum  there,  and  hence  w”(0)  < 0. 

Nov/  the  first  integral  of  Eq.  (6.3)  is 

p 

w"  + 4w  - ew  + C = 0. 


This  equation,  together  v/ith  w(0)  = 1 and  w"(0)  = 0 
gives  C > e - 4. 

Hence,  we  have  for  e = 6: 

a.,  > 1,  a2  < - 1 

and  the  inequality  in  (6.12)  is  again  satisfied. 
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(c)  e > 6,  special  case  e = 120. 

The  procedure  of  cases  (a)  and  (b)  leads  to  a^  > 1 and 
&2  < - 2,  so  that  our  integration  formula  again  holds. 

Since  Eq.  (6.11)  gives  ©0  in  terms  of  e and  C,  the  pro- 
cedure followed  in  the  numerical  computation  is  to  pick  a value 
of  C for  a given  e and  to  find  the  corresponding  ©0.  It  is  due 
to  this  that  the  calculations  for  the  three  values  of  e were  not 
carried  out  for  one  and  the  same  wedge  an<-le,  but  rather  foT 
throe  angles  which  are  only  approximately  equal. 

Wall  Drag.  The  formula  for  the  drag  can  immediately 
be  obtained  from  (6.6)  by  putting  w = 0.  This  gives 

w * ( - ©0)  = - e + 6- . (6.14) 

The  drag  is  readily  obtained  from  this  by  recalling  that  u(r,©)  = 

- Hu  w ( © ) . 
r 

6.2.  The  linearized  result.  In  this  section  we  de- 
note the  radial  velocity  component  by  1 + f(©)  and  Eq.  (6.3) 
becomes 

, f"»  + (4  - 2ce)f » = 0 (6.1 5) 

and  the  boundary  conditions  are  f(0)  = 0 and  f(+  ©0)  = - 1.  We 

9 2 

denote  4 - 2ce  by  and  distinguish  the  three  cases,  s <0, 

o p 

s = 0,  s >0.  The  appropriate  solutions  of  Eq.  (6.1 5)  are 
(a)  s2  > 0, 

cos  s©  - cos  S©0 

W = 1 + f = (6.16) 

1 - COS  S©0 
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(b)  s2  = 0, 

w = 1 ♦ f = 

(c)  s2  <0, 

cosh 

w = 1 

+ f = 

where  s'  = - s. 

The  velocity  gradients 

diately.  They  are: 

(a)  e < -2; 

c 

< 

1 

0 

0 

II 

(b)  e = 

< 

1 

0 

0 

II 

(c)  e >-2; 

< 

1 

<0 

O 

II 

_si 

On 
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(6.17) 


1 - cosh  8 *0, 


(6.18) 


s sin  s0o 
1 - cos  sO, 


3*  sinh  s*0o 
1 - cosh  s'0o 


(6.19) 

(6.20) 
(6.21) 


6.3*  Boundary  layer  solution.  As  a matter  of  interest 
we  consider  the  boundary  layer  solution  of  this  problem.  We 
start  from  Eq,  (6.3)  and  integrate  it  once  to  obtain 


w"  + 4w  - ev^  + C = 0. 


(6,22) 


In  the  exact  solution  we  determined  C from  the  condition 

that  w(0)  =1.  In  the  boundary  layer  solution  we  demand  only 

that  w approach  unity  asymptotically  as  we  go  away  from  the  wall. 

Hence  C can  only  be  determined  in  the  limit  as  e — > 00. 

In  the  interior,  we  expect  w and  its  derivatives  to  be 

of  order  one;  hence,  for  sufficiently  large  e,  Eq.  (6,22)  becomos 

approximately  _ 

- ew^  + C = 0, 


(6,23) 
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N ^«it>  C.  r\r\ 

• • V*«A  -wy  4. 

Since  w is  to  be  one  on  0 -•  0,  we  see  that  C must  be 
equal  to  e within  the  order  of  approximation  of  the  boundary 
layer  theory. 

Following  the  usual  procedure  in  boundary  layer  theory 
we  transform  to  a new  independent  variable 


9 = e1/2(©0  - ©), 


w(9)  = 1(9) 


(6.2'*+) 


in  which  x and  its  derivatives  are  all  of  order  one.  We  have 


ex 


'*  + lfx  - ex  + e = 0. 


(6.25) 


For  e approaching  infinity  this  gives  the  approximate 
boundary  layer  equation 

x"  - x^  + 1 = 0.  (6.26) 

The  boundary  conditions  w(©)  = 1 and  w(©Q)  = 0 transform  to 

lim  x = 1 and  x(0)  = 0.  (6.2  7) 

9 — >00 

One  integration  of  (6.26)  gives 

a = D.  (6.28) 

To  evaluate  D we  apply  the  condition  that  x — >1  as  9 — >co  , and 
hence  x1— ^0  as  9— >cc.  This  gives  D = 2/3.  We  can  now  write 

(x»)2  = | (1  - x)2(x  + 2)  (6.29) 

and 

9 
% 

dn  = 

V 


(6.30) 


x 

n 


0 


dr 


(1  - 0 ^(2  + 0 


r 

- \! 


3 
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Upon  integration  and  substitution  of  w and  © in  terms 
of  't  and  9 this  leads  to  the  following  solution: 

©0  - © = \j&  [tanh"1  \j*  -t-  - tanh'1  \j. 2].  (6. 3D 

The  velocity  gradient  at  the  wall  can  be  found  from 
Eq.  (6*29)  by  putting  'i  - 0.  Hence 

(t»(0))2  = 4/3 

and  „ _ 

w'(-  «„)  = -2-  vT  . (6.32) 

V3 

The  most  important  test  for  applicability  of  the  various 
approximate  theories  is  probably  to  be  found  in  the  values  they 
predict  for  the  velocity  gradient  at  the  wall.  An  inspection  of 
the  table  of  results  shows  that,  for  the  problem  at  hand,  both 
the  standard  and  the  modified  Oseen  technique  give  results  that 
are  accurate  to  within  a few  percent  for  e = 1 and  e - 6,  but 
that  the  modified  method  is  appreciably  better.  For  e of  the 
order  100,  the  Oseen  methods  are  quite  inaccurate,  but  the  bound- 
ary layer  solution  is  already  fairly  good.  It  is  interesting  to 
note  (see  Fig.  6)  that  the  velocity  profile  as  calculated  by 
boundary  layer  theory  is  very  close  to  the  exact  profile  in  a 
region  a tax' ting  at  the  wall  and  extending  to  nearly  one  third  the 
total  channel  width. 

7*  Conclusions.  It  is  evident  that,  in  the  foregoing 
moderate  Reynolds  number  problems,  a successful  prediction  of  the 
macroscopic  feature  of  the  flow  are  given  by  the  linearized 
theory  with  c = .43.  It  can  safely  be  anticipated  that  an 
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equally  successful  application  can  be  anticipated  for  other  flows 
of  a similar  nature.  For  more  complicated  problems,  e. g. , the 
flow  past  a plate  of  finite  chord  at  non-^ero  incidence,  no  con- 
clusion can  be  drawn  until  such  flows  have  been  investigated. 
Similarly,  it  is  not  clear  whether  the  range  of  applicability  can 
be  increased  by  finding  a c(e)  for  larger  e than  those  considered 
here. 


(3)  is: 
and 


8.  Appendix.  The  boundary  value  problem  of  section 
AA4>  - cA^  = 0;  (2.9) 

t|>(x,0)  = 0;  ij>  (x,0)  = - 1 for  x > 0. 


The  solution  as  given  in  [3]  is  found  by  using  Fourier 
transform.  We  define 


00 

n 


-00 


e-i(ot>V)+(Jt>y)dxdy> 


(8.1) 


and  note  that  ^ and  'I’yyy  are  continuous  across  the  plate. 
We  denote  by  f(x)  the  Jumps  in  ^(XjO).  Using  the  foregoing, 
Eq.  (2.9)  becomes 


[(^  + rfV  * ic«^  + q2)]?  = - irf?CO.  (8.2) 


Equation  (8.2)  may  conveniently  be  thought  of  as  the  limit  as 
k— ^0  of  the  equation 

( K 2 + q2  + k2)(q2  + U + ilU  - iklH  = - iqf.  (8.3) 


We  may  also  note  that  the  boundary  condition  on  t^y  at  the  plate 
is  the  limit  as  a— of  \1>  (x,0)  = e“ax  for  x > 0.  We  may  now 
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define 


?U,y)  = 4>(^,n)elny  dy 


2n  I 
u 

-CD 


and  obtain 


<p(S,y)  = 


|y|f[exp{,-  |y|U2+k2)1/Sj-axp{-  | y|  (^l)1/2(c-H01/2j] 


y 21(1  - k)(£  - ik) 


(8.4) 


We  may  now  define  <p^(C,0)  = uQ(x)  x u^(x)  where  uQ(x)  = emax 
where  x > 0 and  uQ(x)  = 0 when  x < 0. 

We  may  now  use  the  usual  Wiener-Kopf  arguments  to  find 
f(0  [the  details  a-e  to  be  found  in  [3]  1 and  find,  when  k— >0, 

a — - 1/2 

f = 2iUA)  * (8.5) 

We  now  can  invert  Eq.  (8.4)  using  this  form  for  f and 

we  obtain 

— -1/2  2 

\|>(x,y)  + y = 2£[tj  erf  (n  Vc)  - (*c)  (1  - e**cT1  ) ] 

= ^g(n)*  (8.6) 

1/2 

Here,  £ + it)  = (x  + iy)  and  \}>  + y is  the  stream  function  for 
the  complete  flow.  The  bracket  is  denoted  by  g(q)  for  conve- 
nience^. 

Consistent  now  with  the  motivation  of  section  (1)  we 
would  like  to  choose  c so  that  this  solution  agrees  with  the  solu- 
tion of  the  precise  problem  both  near  the  plate  leading  edge  and 
far  from  the  plate,  Our  only  information  is  that  for  large  dis- 
tances downstream  of  the  leading  edge  (x  > 20,  say),  the  far 

^ This  explicit  result  (i.e. , 


Eq.  (8.6))  was  not  included  in  [3] 
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field  solution  is  given  by 

£ f(t))  (8.7) 

where  f(r])^2r)  - 1.72  + ...  • The  remaining  terms  do  not  corres- 
pond to  a harmonic  function.  The  asymptotic  behavior  of  ^ = Sg(n) 
becomes  the  same  as  that  of  Eq»  (8.7)  provided  we  choose  c = .4-3» 
Thus,  if  any  function  c(0  is  to  accomplish  the  purpose  outlined 
in  section  (1),  its  value  as  e ~>oo  should  be  .4-3*  In  Fig.  (8.1) 
the  function  f'(n)  and  g r ( n ) are  plotted.  One  should  note  that 
If'/g'  - l|  is  never  greater  than  .1,  the  most  serious  discrepan- 
cy occurring  at  the  plate  surface. 


exact 

solution 

Standard 

Oseen 

Modified 
Oseen 
c = .43 

Boundary 

Layer 

6.63 

6.6/ 

s / A 

O.Otf 

— 

6.00 

6. 26 

5.88 

— 

13.11 

15. 79 

10. 47 

12.65 

Nonr-6 53-00/1 


25 


Bibliography 


[1]  G.  G.  Stokes,  On  the  Theory  of  the  Internal  Friction  df 

Fluids  in  Motion*  Trans.  Cambridge  fhil.  Soc.  , vcl.  8 
pp.  287-319  (W?). 

[2]  Handbuch  der  Ebcperlmentalphyafk.  vol.  4,  part  1,  pp.  2 24- 

31^  (Leipzig,  19315. 

[3  ] J.  A.  Lewis  and  G.  F.  Carrier.  Some  Remarks  on  the  Flat 

Plate  Boundary- layer.  Quart.  Appl.  Math. , vol.  7,  pp. 
228-234  (1949;. 

[4]  G.  F.  Carrier  and  C.  C.  Lin,  On  the  Nature  of  the  Boundary 

Layer  Near  the  Leading  Edge  of  a Flat  Plate,  Quart, 
of  Appl.  Math.  vol.  Vi,  No.  1,  April,  19^8. 

[5]  h.  S.  Ling,  Translational  Motion  of  a Rigid  Circular  Cvllnde 

in  a Viscous  Fluid,  A thesis  submitted  in  partial 
fulfillment  of  the  requirements  for  the  Degree  of 
Master  of  Science  in  the  Graduate  Division  of  Applied 
Mathematics  at  Brown  University,  September,  1949. 

[6]  N.  B.  Haaser.  The  Viscous  Flow  Past  a Flat  Plate,  A thesis 

submitted  in  partial  fulfillment  of  the  requirements 
for  the  Degree  of  Doctor  of  Philosophy  in  the  Graduate 
Division  of  Applied  Mathematics  at  Brown  Univers.ty, 
August,  1950. 

[7]  Private  communication  --  details  expected  to  appear  as 

thesis  at  Brown  University. 

[8]  L.  Bair  stow.  B.  M.  Cave,  and  E.  D.  Lang,  The  Resistance  of 

a Cylinder .Moving  In  a Viscous  Fluid.  Phil.  Trans. 

Roy.  Soc.  London,  vol.  223A,  pp.  383-432  (1923). 


fttyrMld*  number  Ft • 


FIG.  5.2:  Velocity  gradient  ot  tht  plots  for  o*4.  The  dotted  curve 

represents  the  linearized  result.  The  solid  curve  is  the 
finite  difference  ^n- linear  result. 
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